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Abstract 

The phase space structure of certain quantum states reveals structure on a scale that is small compared to 
the Planck area. Using an analog between the wavefunction of a single photon and the electric field of a 
classical ultrashort optical pulse we show that spectral shearing interferometry enables measurement of such 
structures directly. Thereby extending the idea of Praxmeyer et al. In particular, we use multiple-shear 
spectral interferometry to fully characterize a pulse consisting of two sub-pulses which are temporally and 
spectrally disjoint, without a relative-phase ambiguity. This enables us to compute the Wigner distribution 
of the pulse. This spectrographic representation of the pulse field features fringes that are tilted with respect 
to both the time- and frequency axes, showing that in general the shortest sub-Planck distances may not 
be in the directions of the canonical (and easily experimentally accessible) directions. Further, independent 
of this orientation, evidence of the sub-Planck scale of the structure maybe extracted directly from the 
measured signal. 



1. Classical and quantum wave fields 

The electromagnetic field is a proxy for the wave- 
function of a single photon [ij. That is, a spatio- 
temporally localized optical pulse field plays the 
role of a probability density function (with suit- 
able normalizatioii) for the measurement of a single 
quantum of light [2] . Interference phenomena that 
are observed with classical fields therefore bear a 
close relationship to quantum interference phenom- 
ena for single particles, and this provides a route to 
exploring quantum phase space structures in a sim- 
ple way. It is not perhaps surprising that there is a 
similarly close relationship between the characteri- 
zation of a classical ultrafast optical pulse (i.e. the 
measurement of its electric field) and the quantum 
state tomography of a single photon, and this ana- 
log may be used to apply recently developed mea- 
surement techniques in ultrafast optics to the study 
of quantum phase space. Krzysztof Wodkiewicz 
exploited this analog to effect in a number of di- 
rections, most recently in a study of so-called sub- 
Planck phase-space structures discovered by Zurek. 

The key elements of state tomography for lo- 
calized photons are the same as those required to 
fully characterize the electric field of a short optical 
pulse. The central idea for this is related to the 
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Wodkiewicz-Eberly notion of the "physical spec- 
trum" , which is simply the time-resolved intensity 
of a spectrally filtered pulse Q , 



dt'K{t-t';Lu)E{t') 



(1) 



with E{t) the input pulsed electric field, and K{t;ijj) 
the time-stationary causal response function of the 
linear filter. Here, uj is the passband frequency of 
the filter, which acts as a spectrometer, and can 
be tuned. This formulation provides a basis for de- 
veloping methods for measuring E{t), even using 
slow detectors ^j. In fact, the signal, considered as 
a function of the two variables t and w is a sono- 
gram of the input pulse. That is, a positive-definite 
phase-space representation of a field (as we consider 
only pulsed sources in this paper, we shall assume 
the input field has compact support), obtained by 
selecting a narrow spectral component of the pulse, 
and measuring its time-dependent intensity after 
the filter. By adjusting the passband frequency and 
measuring the temporal intensity at each frequency, 
a time-frequency representation - the sonogram - is 
built up. 

More generally, the fast-responding detector that 
is required for measuring Eberly- Wodkiewicz phys- 
ical spectrum can be replaced by a slow, time- 
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integrating detector in front of which there is a 
time-non-stationary Unear fiher, such as a fast shut- 
ter. This may be modehed as a hnear filter with 
response function G{t; r), such that the fiher trans- 
mission is unity for a smaU time window near t = t. 
The transmitted field after this "time gate" is 



Eoutit) = G{t;T)E,„{t). 



(2) 



A more common approach to the characterization 
of ultrashort optical pulses involves these two ele- 
ments in reverse order. That is, a temporal slice of 
the test pulse field is taken using a shutter, and the 
transmitted spectrum measured using a spectrom- 
eter and a slow detector. The resulting positive 
definite time-frequency representation of the pulse 
is a spectrogram. In fact, if the spectrometer pass- 
band is very narrow, the detected signal is closely 
related to a Gabor transform 



dtG{t - T)E{t)e- 



(3) 



By adjusting the passband frequency and the delay 
T, over a sufficient range, it can be shown that the 
detected signal is tomographically complete. That 
is, it can be inverted to extract E{t). It is easy 
to see from the form of the signal that it has the 
structure of a phase-space representation. Indeed, 
if G{t) is a Gaussian, the functional form is just that 
of the Husimi distribution in quantum physics (or 
the Q-function in quantum optics) for a pure quan- 
tum state. Therefore one might expect interference 
effects parallel to those of interest in quantum sys- 
tems to appear in the spectrogram. 

One approach to measuring the electric field 
based on spectrograms is frequency-resolved opti- 
cal gating (FROG) [|, in which the "gate function" 
of ([3]), G(t), is implemented by means of a nonlin- 
ear optical interaction. The test pulse field E{t) 
is commonly mixed with a time-delayed replica of 
itself in a material with a x^^^ nonlinear response. 
This generates the nonlinear spectrogram 



5(t,^) 



dtE{t - T)E{t)e'' 



(4) 



Extracting the field from this nonlinear spectro- 
gram requires an iterative deconvolution algorithm 
[6|. However, for the purposes of the measurements 
required in what follows this is not necessary. 

The Gabor spectrogram is a particular example 
of a general scheme for pulse characterization that 



involves a sequence of two linear filters. In that 
case, it is a time-non-stationary amplitude filter (ef- 
fectively a shutter) followed by a time-stationary 
amplitude filter (a spectrometer), as depicted in 
Fig. [IJa). There are four possible minimal se- 
quences comprised of amplitude and phase-only fil- 
ters. These may be categorized both by the par- 
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Figure 1: Linear filter representation of two pulse character- 
ization schemes, (a) Gabor spectrogram, consisting of time- 
non-stationary gate G{t; r) followed by spectral filter H{lu). 
(b) Spectral shearing interferometry, with a phase modula- 
tor producing a frequency shift P(t; Q), delay line D(t) and 
spectral filter H(uj). 

ticular configuration of filters and by the algorithm 
used to extract the field from the signal ^] . 

A complementary approach is to use the filters in 
parallel. There are similarly four possible configu- 
rations 3 of apparatus in which the input pulses 
are split into two at a beamsplitter, with each com- 
ponent traversing a filter, one of which is time- 
stationary, and the other time non-stationary, after 
which the filtered pulses are recombined and the 
spectrum of the interfering pair registered. A par- 
ticular example of this category is spectral shear- 
ing interferometry (SSI). In this case one of the 
filters following the beamsplitter is a phase mod- 
ulator. The role of this non-stationary phase-only 
filter is to shift the frequency of the pulse by an 
amount H. — the spectral shear. The parallel filter 
has a time-stationary response, and causes a shift 
in the arrival time of the second pulse replica — it 
is a delay line. The pulses are recombined on a sec- 
ond beamsplitter and the joint spectrum measured. 
The linear filter representation of SSI is given in 
Fig. [ijb) . The resulting spectral interference pat- 
tern is 



S{uj; T, n) = \E{uj) + E{lo - r2)e- 



(5) 



Again, when the modulator inducing the frequency 
shift is synthesized using nonlinear optics, this ap- 
proach is called spectral phase interferometry for 



direct electric field reconstruction (SPIDER) 
The field reconstruction in this case is linear and 
non-iterative. It is based on etracting the spec- 
tral phase function of the field (f){uj) = ArgE{uj). 
A simple Fourier transform of the signal S{oj; r, fl) 
with respect to ft enables the interference terms 
E*{ui)E{uj — ri)e~'"'^ to be separated, and the ar- 
gument of this term, with a reference phase sub- 
tracted, is simply T[ijj) = (/)(cij) — — fJ). This 
sampling of the phase function is sufhcient to en- 
able the pulse field to be accurately estimated on 
the temporal domain [— Tr/fJ, 7r/il]. 

Using both FROG and SPIDER is it is possible 
to construct a Wigner phase space representation 
of the pulse directly. It is this function that is com- 
monly used in quantum mechanics to study non- 
classical properties of quantum states. Of particu- 
lar interest are the types of "sub-Planck" structures 
discussed by Zurek [lO| . He showed that two states 
with significant overlap in phase space may be in 
fact orthogonal, even if they are separated in either 
dimension by a degree that is significantly smaller 
than the scale implied by Planck's constant. 

Using the measurement methods discussed 
above, it is also possible to extract directly a signal 
component that illustrates the important structure 
of the field. This is a considerably more efficient ap- 
proach than a full characterization of the field, or, 
equivalently for single photons, full state tomogra- 
phy. In effect, these components are "witnesses" to 
the sub-Planck structure of the field. In the next 
section we introduce sub-Planck structure and de- 
scribe the measurements necessary for implement- 
ing witnesses. 



2. Sub-Planck structure 

One measure of the sensitivity of a quantum state 
to perturbations is the size of the smallest displace- 
ment — {5x, Sp) in phase space required to make 
the state orthogonal to the original; that is, for the 
inner product 



(V' I i''} 



between state | i/;) and its displaced replica | i/;') 
to be zero. Zurek showed that this minimum dis- 
placement along a given direction in phase space 
is inversely proportional to the spread of the state 
in the perpendicular direction. For example, the 



lowest-order estimate of the "zero-overlap" momen- 
tum shift 6ri is 



(7) 



i.e. inversely proportional to the root-mean square 
spread of the state. The size of the zero-overlap 
shift can be significantly smaller than the Planck 
scale y/h. This is connected to the existence of 
phase space structures of dimension given by ([7]) 
which appear in the Wigner representation of the 
state, 



1 

27rn 



) z z 

(8) 

By means of arbitrary rotations of phase space, 
this argument can be applied along any direction 
in phase space: the smallest displacement leading 
to orthogonality and the sub-Planck structure both 
have a size inversely proportional to the spread of 
the wave function along the perpendicular direc- 
tion. 

Some examples of Wigner distributions are given 
in Fig. [21 We use atomic units, with h — 1. 
Fig. Ufa) shows a wavefunction consisting of a su- 
perposition of two Gaussian states separated posi- 
tion by xo = 4. Sub-Planck structure is evident, 
in the form of interference fringes along the mo- 
mentum axis. The magnitude of the overlap as a 
function of (Sx, 6p) is shown in Fig.[2Jb). The over- 
lap also contains fringes, modulated proportional 
to 1 + cos{Spro). The smallest shift which produces 
a zero-overlap, and hence orthogonality, occurs for 
Sp — ±7r/ro. For the clarity of Fig. [21 we did not 
choose ro to be particuarly large, but it is clear that 
as ro increases both the fringe period and area can 
decrease well below their respective Planck length 
and area scales. 

The sub-Planck structure of the the analogous 
classical optical field, with time and frequency sub- 
stituted for position and momentum, has been stud- 
ied by Wodkiewicz and coworkers ^11] . The Wigner 
distribution of a double pulse, consisting of two 
identical sub-pulses with time delay Iq, contains 
spectral interference fringes of period 27r/to To di- 
rectly observe the scalar product, Wodkiewicz and 
others pointed out that the spectrum of the second 
harmonic of an ultrashort pulse 
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Figure 2: Wigner and overlap functions of quantum states 
consisting of two separate Gaussians. (a) Wigner function, 
Gaussian states separate in position. The blue lines show 
the marginals, (b) Overlap function corresponding to (a), 
(c) Wigner function with the Gaussian states separate in po- 
sition and momentum, (d) Overlap function corresponding 
to (b). 



is equivalent to the scalar product of the pulse with 
a frequency shifted replica 



\{E{t)\E{t), 



E*{t)E{t)e 



'.Scot 



(10) 



when the electric field has flat temporal phase so 
that the presence of the complex conjugate in pO|) 
makes no difference. Using Frequency-Resolved 
Optical Gating (FROG) they directly observed the 
scalar product for various frequency displacements. 

This arrangement enables the direct observation 
of the overlap between a pulse and its frequency- 
shifted replica, with the additional proviso that the 
pulse has a linear temporal phase. In many cir- 
cumstances it is desireable to measure the overlap 
with arbitrary linear combinations of time- and fre- 
quency shifts, and for arbitrary pulses. FiglDJc) 
shows a superposition of Gaussian states which 
are displaced in both position and momentum by 
To = {xo,Po) = (4,8). Notably, although the indi- 
vidual states do not overlap in either domain, in- 
terference fringes appear in the Wigner distribution 
aligned parallel to the separation vector. The am- 
plitude of the overlap is shown in Fig. [2^c). In this 



domain, the central fringe modulation is propor- 
tional to 1 + cos(Jq • ko), where the wavenumber 
ko = {j>Oi~xo) is perpendicular to the separation 
Yq. The smallest shift which produces a zero is 
±7rko/|kop = ±7r(-po,a;o)/(a;o + pI), again per- 
pendicular to Vq. 

We shall now show that the scalar product be- 
tween a pulse and a time- and frequency-shifted 
replica is given by the same raw data which is used 
in SSI. As described in Section [T] above, in SSI the 
interferometric term 



D{uj,n) ^ E*{uj)E[uj-n) 



(11) 



is obtained via Fourier filtering of the interferogram 
and removal of the carrier phase. For normal pulse 
reconstruction, the phase of D{lo, 51) is used to re- 
construct the spectral phase of the unknown pulse. 
However, here we note that by setting f2 = (5^, the 
Fourier transform of D{uj,5t^) along lo is given by 

/oc 
E{t)E*{t-5t)e-''-'dt (12) 
-OO 

= {E{t-5t)e''-'\E[t)) (13) 

which is exactly the overlap between the field and 
replica shifted in frequency and 5t in time. 
Therefore, the Fourier transform of the filtered side- 
band provides the overlap directly without any in- 
termediate processing. Zeros of D(dt,S^) around 
the origin are witnesses to the sub-Planck structure 
of the quantum state analogous to the unknown 
pulse. 

In the next section, we describe an experimen- 
tal setup for observing witnesses to the sub-Planck 
structure in an ultrashort pulse analagous to the 
position- and momentum-separated quantum state 
in Fig.djc). 

3. Experiment 

There are many implementations of SSI that can 
be used to measure D(uj,di^) for single frequency 
shifts. For convenient, single-shot measurement of a 
continuous ran ge o f frequency shifts, we used SEA- 
CAR-SPIDER [13], a recent twist on SSI that per- 
mits the simultaneous acquisition of multiple shears 
and a precise and robust calibration procedure. The 
principle of SEA-CAR-SPIDER is shown in Fig. H 
Two spatially chirped fields, known as the ancillae, 
are prepared and sum-frequency mixed with the un- 
known pulse in a crystal. The the ancilla beams 




Figure 3: SEA-CAR-SPIDER concept; the unknown pulse 
(UP) undergoes sum-frequency mixing with two oppositely 
spatially chirped ancillae (Ai, A2) in a x'^' crystal. The 
upconverted beams are re-imaged onto a two-dimensional 
imaging spectrometer. 



are oppositely spatially chirped, so that in the crys- 
tal their local frequencies at transverse position x 
(indicated in Fig. ^ are Wup ± olx. Here, cjup is 
their common upconversion frequency at x = and 
a is the amount of spatial chirp. Two upconverted 
replicas of the unknown pulse are produced, with 
a relative spectral shear of S^^ = 2ax. These are 
re-imaged onto an imaging spectrometer, their con- 
verging angle producing a spatial interference pat- 
tern 

S{ijj + u!^^p,x) = \E{ijj + ax)e^^'''^ + E(u} - ax)\^ 

(14) 

where kc is the wavenumber of the interference 
fringes, determined by the beam convergence an- 
gle. To process this interferogram, either for re- 
construction of the phase of the unknown pulse or 
for determination of the overlap (|13p . we apply the 
two-dimensional Fourier transform to the (i, kx) do- 
main. Here, the transforms of the spectral intensi- 
ties \E{bj ± aa;)^ appear at the origin, whilst the 
interference term between the replicas manifests as 
two sidebands located at (0, ±fcc). One of these is 
isolated and the inverse transform applied. Remov- 
ing the spatial carrier kcX and using the linear rela- 
tion between x and the shear, we obtain D{uj,5ui). 

The optical pulse analogue to the quantum state 
of Fig. [5] is a bichromatic double pulse — i.e. 
two sub-pulses with both a temporal and spec- 
tral separation. We synthesized these using a 4-f 
pulse shaper with a « 1.1mm thick glass slide 
in one half of the Fourier plane, so that frequen- 
cies above 2.287radfs~^ were delayed. We also 
placed a thin opaque block at the edge of the 
glass slide, attenuating frequencies from 2.284 — 



— 2.291 rad fs^^. Through this procedure, we cre- 
ated a bichromatic double pulse with temporal sep- 
aration to = 1820 fs and spectral separation luq = 
18mradfs~^, which we then characterized using 
the SEA-CAR-SPIDER apparatus. 

The raw data trace is shown in Fig.lH The fringe 
pattern can be interpreted as follows: the phase of 
each fringe is modulated by (I)(lu + S^j) — 4>{ijJ — 5^^)- 
To either side of the spectral gap at 2.287 rad fs~^, 
the two frequencies uo both reside in the same 
sub-pulse. Each sub-pulse is transform limited, and 
hence the phase difference between the two frequen- 
cies is a constant. The fringes are therefore unmod- 
ulated away from the spectral gap. Closer to the 
spectral gap, and for larger shears, m ± 5^ reside 
in different sub-pulses. Their different group delay 
produces a phase difference which varies linearly 
with frequency. The fringes are therefore tilted in 
this region. 
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Figure 4: Raw SEA-CAR-SPIDER data for the bichromatic 
double pulse. 



Although our aim was to directly measure the 
overlap (fT^ . we also reconstructed the pulse to 
verify our synthesis and to numerically calculate 
its Wigner function. Because the pulse contained 
significant gaps in its spectrum, it was necessary 
to use the multiple shear reconstruction algorithm, 
[ist . a recently developed generalization of stan- 
dard SSI concatenation procedure Q which en- 
ables the use of a large spectral shear to bridge 
spectral gaps. The measured spectrum and tem- 



poral profile of the bichromatic double pulse is 
shown in Fig. \E[ We calculated the Wigner dis- 
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Figure 5: (a) Spectral amplitude (solid, left axis) and phase 
(dashed, right axis) of the bichromatic double pulse, (b) 
Corresponding temporal intensity (solid, left axis) and phase 
(dashed, right axis). 

tribution of the measured pulse, shown in Fig. [HI 
The two peaks at {t,uj) = (-423 fs, 2.275 rad/fs) 
and (1387 fs, 2.293 rad/fs) represent the sub-pulses, 
whilst the interference fringes between them are 
aligned parallel to their displacement (to,t<^o) in 
phase space. As expected, Fig.|n]bears a close qual- 
itative similarity to Fig.[2^c). 

We calculated the overlap of the pulse and its 
phase-space-displaced replica D{5t,5S} by Fourier 
transforming D{u!,di^). The result is shown in 
Fig. [3 and as expected bears a close similarity 
to Fig. E^d), with fringe modulation proportional 
to 1 -I- cos [{6uj,St) ■ (wo, — ^o)]- Zero overlap occurs 
whenever {Suj,5t) ■ (wo, -^o) = iT^. 



4. Summary and conclusions 

The analogy between the electric field of an ul- 
trashort pulse and the wavefunction of a single pho- 
ton can be exploited so that classical wave inter- 
ference in optics mimics quantum interference, en- 
abling observation of phase-space structures that 
have important implications for dissipative quan- 
tum dynamics. 



Sub-Planck structure is just such a quantum in- 
terference effect which determines the sensitivity of 
a quantum state to perturbations, considered as dis- 
placements in phase space. We have demonstrated 
that spectral shearing interferometry directly mea- 
sures the overlap between an ultrashort pulse and a 
replica that has displaced in phase space in a range 
of directions, thereby providing a direct test of the 
Zurek hypothesis. 

An important question is whether it might be 
possible to implement the same approach for a 
quantum system consisting of massive particles? In 
fact, it is possible to devise a shearing interferome- 
ter for the center-of-mass wavefunction of an atom. 
[3] A momentum shift of the transverse degree of 
freedom can be accomplished by passing the atom 
through a standing-wave light field in the Lamb- 
Dicke limit. Propagation in free space leads to 
the detection of an atom interferogram with sim- 
ilar properties to those of the SPIDER method dis- 
cussed here. 

Such methods may enable a deeper exploration of 
uniquely quantum characteristics, those which dis- 
tinguish purely quantum dynamics from that possi- 
ble in classical systems, to which question Kryszstof 
Wodiewicz devoted a great deal of his work. 
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Figure 7: Phase space autocorrelation of the bichromatic 
double pulse. 
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